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Abstract 

The  method  developed  by  Wiener  and  Hopf  for  solving  a  certain  class 
of  integral  equations  is  applied  to  the  problem  of  a  discontinuity  in 
cross-section  of  a  rectangular  waveguide o  We  consider  the  two-dimensional 
acoustic  problem  of  two  channels  of  infinite  width  and  different  heights 
joined  together  with  a  short  step;  it  is  assumed  that  the  fundamental  mod« 
is  incident  in  the  larger  channel  traveling  towards  the  step.   For  values 
of  the  channel  height  such  that  only  the  lUndemental  mode  i  s  propagated 
without  attenuation,  this  problem  has  been  treated  in  the  Wave  Guide  Hand- 
book; our  results  are  in  very  good  agreement  with  the  results  given  there. 
The  present  analysis  is  extended  to  cover  the  case  where  the  larger  channel 
is  capable  of  supporting  the  mode  above  the  fundamental  without  attenuationi 
Some  numeirical  values  are  given  for  the  ratio  between  the  reflected  and  in- 
cident energies  for  various  values  of  the  channel  and  step  heights. 


Table  of  Contents 


Page 


I.  Introduction  1 

II.  Statement  of  the  problem  3 

III.  Solution  for  T  U 

IV,  Calculation  of  reflection  coefficient  10 

V,  Comparison  with  the  results  of  the  Handbook  13 

VI.  The  case  of  ka/n  >  1  l5 

VII,  Solution  of  the  equations  16 

VIII,  Numerical  results  23 

Appendix  A  26 

Appendix  B  28 

References  32 


-  1  - 

I.   Introduction 

Problems  involving  discontinuities  in  waveguides  are  of  great  practical 
importance  and  a  considerable  amount  of  attention  has  been  devoted  to  them  in 
recent  3rears,  The  two  types  of  discontinuities  that  have  been  most  extensive- 
ly treated  for  a  rectangulair  guide  are  the  bifurcation  of  the  guide  by  a  semi- 
infinite  plane  and  a  step  discontinuity  in  the  guide  cross-section.  The  problem 
of  the  bifurcation  may  be  solved  in  terms  of  known  functions  by  applying  the 
well-known  Wiener-Hopf  technique  for  solving  integral  equations.  The  methods 
for  solving  the  problem  of  a  step  discontinuity  which  have  appeared  in  the 
literature  are,  however,  in  no  way  based  on  the  Wiener-Hopf  approach.   Two 
distinct  methods  have  been  employedj  the  first  consists  of  expanding  the  field 
components  on  either  side  of  the  discontinuity  in  terms  of  the  appropriate 
eigenfunctions  for  the  region  considered,  and  then  matching  the  components 
across  the  plane  of  discontinuity.  The  second  method  used  is  the  equivalent 
static  approach  as  described  in  the  Wave  Guide  Handbook'-  •' ,  Since  there  is 
obviously  a  physical  relationship  between  these  tvjo  problems  it  is  desirable 
to  solve  them  both  by  the  same  method  in  order  that  the  solution  should  dis- 
play clearly  the  relationship  between  them,  and  this  is  the  object  of  the  pre- 
sent work. 

We  therefore  consider  the  application  of  the  Wiener-Hopf  technique  to 
solve  the  problem  of  a  step  discontinuity  in  a  parallel  plate  region  (or, 
equivalently,  the  coupling  of  two  parallel  plate  waveguides  of  different 
height).  The  problem  is  formulated  as  an  acoustic  one,  but  it  could  have 
been  equally-well  stated  in  electromagnetic  terms  by  considering  the  case  in 
which  the  only  non-vanishing  component  of  the  magnetic  field  is  parallel  to 
the  edges  of  the  step.  It  is  assumed  that  the  fundamental  mode  is  incident 
from  the  larger  region  on  the  step  and,  initally,  that  the  dimensions 
of  both  regions  concerned  are  taken  to  be  such  that  only  the  fundamental 
mode  may  propagate  without  attenuation.  This  restriction  is  later  removed  and 


-  2  - 


the  case  of  the  larger  region  being  capable  of  sustaining  two  non-attenuated 
modes  is  treated.  The  advantage  of  this  method  of  attack  is  that  it  affords 
a  clearer  physical  pictui^  of  the  structure  of  the  solution  than  that  obtained 
by  the  other  methods.  In  particular,  it  is  possible  to  separate  the  reflection 
coefficient  for  the  fundamental  mode  into  two  parts,  cne  of  which  is  expressible 
in  terms  of  known  functions  and  another  which  may  be  obtained  from  the  solution 
of  an  infinite  set  of  equations.  That  part  of  the  reflection  coefficient  which 
is  known  exacly  is  by  far  the  most  dominant  part  of  the  complete  expression, 
and  thus  the  percentage  error  in  detennining  an  approximate  value  of  the  reflec- 
tion coefficient  is  much  less  than  the  corresponding  error  made  in  the  solution 
of  the  equations.  The  analytic  part  of  the  solution  may  be  obtained  by  simple 
physical  arguments  from  the  known  solution  of  the  bifurcation  problem  and  thus 
the  present  form  of  the  solution  gives  a  distinct  connection  between  the  two 
types  of  problems. 

For  the  case  in  which  only  the  fundamental  mode  is  propagated  without 
attenuation  in  any  of  the  regions  considered,  it  is  possible  to  compare  the 
present  results  with  those  given  in  the  Wave  Guide  Handbook  which  were  obtain- 
ed by  the  equivalent  static  method,  and  the  results  of  both  methods  show  re- 
markably good  agreement.  The  form  of  the  infinite  set  of  equations  obtained 
here  is  such  that  it  is  possible  to  estimate  the  error  made  in  solving  only  a 
finite  niimber  of  them,  and  a  method  is  indicated  by  which  this  error  could  be 
estimated  in  any  particular  case.  The  method  used  here  may  be  easily  extended 
to  solve  the  problem  of  a  bifurcated  waveguide  when  the  material  composing  one 
of  the  bifurcated  regions  differs  in  dielectric  properties  from  that  of  the  re- 
mainder of  the  waveguide.  The  approach  may  also  be  used  to  consider  the  problem 
of  radiation  from  a  dielectric  waveguide,  and  it  is  hoped  that  both  of  these 
problems  will  be  treated  in  future  papers.  The  method  can, of  course,  be  applied 
to  the  general  electromagnetic  case  since  it  is  known  that  this  problem  can  be 
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reduced  to  the  solution  of  two  scalar  problems'-  J,  but  no  attempt  is  made 
to  consider  this  more   general  case.     The  author  has  been  informed  by  Prof, 
N.  Marcuvitz  that  Prof,  J,  Schwinger  has  considered  the  application  of  the 
Wiener-Hopf  technique  to  problems   similar  to  the  one  treated  here.     An  account 
of  his  work  may  appear  in  a  forthcoming  book  by  Marcuvitz  and  Schwinger, 

The  author  would  like  to  acknowledge  the  permission  of  Mr,  R.A,  Brooker 
to  use  the  digital  computer  at  Manchester  University  for  most  of  the  computa- 
tions involved  in  the  text,   and  also  the  many  helpful  suggestions  given  by 
Mr.   D,   S.  Jones  during  the  period  in  which  the  present  work  was  carried  out. 


II.     Statement  of  the  problem 


d  =  Q-c 


O/////////////////////^  X 


Fig.  I 


The  waveguide  is  assumed  to  occupy  the  region  -oo  <  x  <  0,  d  <  y  <  a;  and 
0  <  X  <  00,  0  <  y  <  a  (see  Fig.  I),  where  x,  y  are  Cartesian  coordinates  with 
their  origin  at  0,  Oz,  the  third  axis  of  the  system,  is  not  shown,  and  the 
solution  is  assumed  to  be  independent  of  z,  A  time  variation  e    will  be 


taken  as  imderstood  throughout  the  work,  and  it  is  assujned  that  an  incident 

wave  e         is  appi-oaching  the  step  from  the  right,  k  being  2  ii/(  wave  length). 

ikx 
The  total  velocity  potential  will  be  written  as  e         +  ▼»  and  hence   a 

function  T  is  sought  which  satisfies  the  following  conditions: 

(i)  (7^  +  k^)Y  -  0. 

(ii)  |I    -  -ik    on    X  -  0,        0  <  y  <  d 

ox  ~      ~ 

(iii)  |I    -     0       on    y  -  a,       -oo  <  x  <   oo, 

on     y  ■  d,       -00  <  X  <  0,   and 

ony"0,  0<x<oo  , 

It  is  convenient  to  assume  that  k  «  k     -  ik.,  where  k.,  k^  >  0;  then  the  form 
of  the  radiation  condition  for  T  is  T  -  0(6""^    '*')   as    |x|  ->oo.     In  order  to 
specify  the  problem  uniquely,  a  condition  has  to  be   imposed  on  T  at  the  edges 
of  the  step  t-  J.     This  condition  is  that  T  is     0(i)  as  the  point  of  observa- 
tion approaches   (0,d)  or  (0,0),  while    \Vf\   is     0(r"  '^)  as  the  point  of  ob- 
servation approaches   (C,d)   and  is     0(1)  in  the  neighborhood  of  the  origin - 
-  r  being  the  distance   from  the  point  approached.     The  above  conditions  now 
specify  the  problem  completely  and  uniquely. 


III,   Solution  for  Y 

We  shall  now  define  a  function  {2)(s,y)  by 

0(s,y)     -  e'®^  T(x,y)dx 

■'-00 


where  s  is  the  complex  variable  (T  +  if.  From  the  assumption  on  T  as  |x|  ->oo, 
it  is  clear  that  the  above  integral  is  convergent  for  -k.  <  0'<  k^^.  Hence  ^  is 
a  regular  function  of  s  in  the  region  |  (T  |  <  k.,  Ve  also  define  5J_(s,y), 
9f+(s,y)  as 
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0  00 

-ioo  "Ti 


Clearly,  from  the  behavior  of  T  as  |x|  ->oo,  ^^(s,y)  is  a  regular  function  of 
s  for  G'<  k . ,  and  ^  (s,y)  is  a  regular  function  of  s  for  (T  >  -k..  In  the 
region  0  <  y  <  d^  9'_(s,y)  is  identically  zero  and  hence  we  have  the  following 
equations.  For  y  >  d 

(1)  ^  +  ^<^     '     0. 

dy 

For  0  <  y  <  d 

2 

(2)  11+  +   k\  -  -ik  +  sf(y) 
dy 

2  2  1/2 
whei^  f(y)  »  T(0,y)  and  K  »  (s  +k  )  '  ,  and  where  that  branch  of  the  root  is 

taken  which  reduces  to  k  when  s  ■  0» 

The  solution  of  Eq,  (2)  which  satisfies  the  boundary  condition  on  y  ■  0  is 


(s,y)  =  A  cos  Ky  -  -- 


sin  K(y-t)dt  +  I 


/ 


f(t)sin  Wy-t)dt 


(3) 


y 
A  cos  Ky  -  ^  (l-cosKy)  +  |  f  f(t)sin  K(y-t)dt. 
K  Jo 


where  A  is  a  function  of  s  but  not  of  y. 

If  primes  are  taken  to  mean  derivatives  with  respect  to  y,  then  differ- 
entiation of  Eq.  (3)  and  then  substitution  of  y  -  d  in  the  resulting  equation 

gives 

d 

(U) 

'0 


(f\{s,d)     «  -  K  A  sin  Kd  -  i^  sin  Kd  ♦  s   f(t)cos  K(d-t)dt. 


If  the  value  of  A  from  Eq.  (U)  is  substituted  in  Eq.  (3),  then  this  latter 
equation  may  b  e  rewritten  as 
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KsinKd  (2^(s,y)-  -  ^^^(6,d)+  ~  sin  Kd-s    (  f  (t)cos  K(c^-t)dt  ^  cos  Ky 


(5) 


*  <  Z  (  f(t)sin  K(y-t)dt  -  ifc  (1-cos  H,y)l  Ksinl^d   . 


By  definition  (jl  (s,y)  is  regular  for    <r'>  k.  «     Hence,   since    K.sin  K>d  has  zeros 

in    Cr'>  -k , ,  the  right-hand  side  of  Fq,    (3)  must  vanish  for  those   zeros  of 

K  sin   K^d  in  C>  -k . ,     These   zeros  will  be  at  s  ■    K      (m  ■  0,1,2, ,..)  where 

2  2 
.  /m  n         ,  2^   1/2 

m  ^2 

and  is  positive  when  real,   and  positive  imaginary  when  not  real   (for  k  real). 
It  is  also  easy  to  show  that  R   K     ^  ^-t  •     Thus,   in  order  that  the  right  hand 
side  of  Eq,   (5)  vanishes  for  s  »    K    »  the  following  relations  must  hold 

iiJ^(K„,d)  -   (-)"'Kj„    I  f(t)cos  JllJ^     dt  m  /  0 


9(^(ik,d)+ikd        / 


TF 


f(t)dt   . 
0 


If  the  expression  for  f(t)  as  a  Fourier  cosine  series  as  obtained  from  the 
above  expressioneis  substituted  in  Eq.  (5),  then  application  of  the  well- 
known  results  for  term  by  term  integration  of  Fourier  series,  together  with 
some  slight  manipulation,  yields  the  following  equation  for  {i'^(s,y) 


^  mny 

(7)  Ksin  Kd9^(s,y)-  -  5i(^(s,y)cosK  y  +  ^^^gSf^  *  ^  sinWd  ^I  F„  -^ 

n=o    s  -  K 

m 

where 
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We  now  consider  the  region  y  >  d.  The  solution  of  Eq.  (1)  which  satisfies 
the  boundary  condition  on  y  =  a  is 

9f(s,y)  =  B  cos  K(y-a) 

where  B  is  a  function  of  s  only.     Since  SY/Qy  vanishes  on  y  »  d  for  -co  <  x  <  0, 
(^  (s,d)  »  0  and  hence 

-B  Ksin  K(d-a)  -  9f^(s,d). 

Thus  we  obtain  the  following  equation  for  ^(s,y) 

rt/       \  COS  K(y-a)       rtV     j\ 

Hence 

(8)  i2^(s,d)  *  9f_(s,d)  -  -  °°^^^'^"^^     ^/s,d). 

Substitution  for  0  (s,d)  from  Eq.   (7)  into  Eq.   (8)  gives 

fc,\  rt  /„  ^^  sin   Ka  ^* ,     a\       '^  2s  i^     ^''  ^m 

^9)  ^Js»d)  -    Ksin  Kc  sin  Kd  ^J^*^^'  F^IT?  '  d"  2_^  X^  * 

m 


On  writing 


sin  K^a 

K^sin  K,c  sin  Kd 


as 

TIJT 

it  is  seen  that  Eq.   (9)  may  be  rewii.tten  as 

N_(s)  00      (-)\ 

N  (s)9J  (s,d)  -  N^(s)^_^(s,d)-  -^^^r^  -  -g-  N_(s)5;;;  -2 — ^  . 

m=o  8  -  K„ 
m 

The  formulas  for  N  (s),  N_(s)  are  given  in  Appendix  A  but  a  brief  summary  of 
their  essential  properties  will  be   given  here.     N  (s)  is  regular  and  non-zero 
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in  the  region  (r>  -k, ,  N  (s)  is  regular  and  non-zero  for  5~<  k . ,     Also,  as 
|s|  -i>oo,  N  (s)'vs'*    '     and  N_(  s )  ^v^  s  '    • 

It  is  now  convenient  to  rewrite  the  above  equation  for  (}  (s,d)   as 
follows 

N  (s)-N   (-ik)       „^  00     (-)"g         ,    00   (-)"n   (-K  )g_ 

m-o  8     -H  m«o     ^         '*in'' 

m 

(10)  ,  N  (-Ik)       ,    00     (-)'\  (-X  )g 


in=o  m 


Clearly  the  left  hand  side  of  equation  (10)   is  regular  for  <r<  k     while  the 
right  hand  side  is  regular  for  6' >  -k. j  hence  both  sides  are  equal  to  an 
integral  function  in  the  coramon  strip   |  <r  |  <  k..     This  integral  function  may- 
be determined  from  examination  of  the  growth  of  both  sides  of  equation  (10) 
as    |s|  ->oo.     In  order  to  do  this  we  require  the  behavior  of  (?'(s,d)  as 
|s|  ->oo  in    (r>  -k^  and  the  edge  condition  shows  that  ^  (s,d)  '^  s'^^^  as 
|s|  ->oo.     The  integral  function  is  now  eaily  seen  to  be  zero  and  we  hava 

,  N  (-ik)       .    oo     (-)"n  (-X  )g„ 

(11)         »,(.)<(.,d) .  ^^^ .  \  n      -\/^ . 

m«o  "^m 

From  Eq.  (11),  on  substituting  s  »  K  ,  the  following  set  of  equations  is 


obtained  to  determine  the  0  ( k  ,d) 

+       m 


1  1   00     N  i'tj\       a^N  (-ik) 

where  a^  -  (jf^i  \<,^,d)  /N^(-ik). 

For  the  case  in  which  only  one  mode  is  propagated  without  attentuation  in 
both  regions  of  Fig.  I,  it  is  possible  to  reduce  the   above   set  of  equations 
to  a  real  set  as  follows:     we  define  L  (s)   »  N  (s)(6+ik),  L_(s)   »  N__(sy(s-ik). 


From  the  formulas  for  N  (s)  and  N_(s)  it  is  clear  that  L  (  K  )  is  real  for 


n  >  1.  Hence  Eq.  (12)  may  be  rewritten  as 


<-    ^*-"^--&^S^<^V^^>-?&^ 


•^-K^    V.v-ik) 


m 
m 


From  the  asymptotic  behavior  of  ^  (s,d)  and  N_(s)  it  is  seen  that  the  terms 
of  the  second  series  in  Eq.  (13)  are  O (m"    )  as  m  ->oo,  hence  this  series 
is  convergent. 

Putting  n  «  0  in  Eq.  (13),  we  obtain 


(iU) 


m=l 


^  \^J'\?  ^I-J-i^) 


K 


If  the  value  for  the  series 

00 


m»l 


as  obtained  from  Eq.  (lU)  is  substituted  into  Eq.  (13)  then  the  following 
system  of  equations  is  obtained 

[^  I  m«l  "^m^  ^m       ^n' 


n  -  1,2,    •••   . 


Ifefining  a  new  set  of  constants  c     for  n  >  i  by 


2a 
n 


"       1*  a^L^(ik)-  5aL_(-ik) 


it  is  seen  that  these  will  satisfy  the  following;  set  of  equations 

00     c  L  (-KjiK^K-k^) 

(16)  V>^nV^-li;     ^'k    (k      .K    )  -       (-i,2,...). 

m«l         '^a^  *m       "•  n' 
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The  set  of  Eq.  (16)  is  real  only  when  the  fundamental  mode  is  propagated  with- 
out attenuation.  Finally  substitution  for  a  into  Eq,  (lU)  yields  the  following 
equation  for  a  . 


L  (ik)a  -1  +(a  /d)L_(-ik) 
(17)  *  P E -     ig 

-L^(ik)a^-1  +(a^/d)L_(-ik) 


where 


,      00     c  L  (-  K) 
k     t —       ra  -         m 


K     t —       ra  - 

'  ■  ^  £i  "~^ 


m 


IV»  Calculation  of  reflection  coefficient 

In  problems  involving  discontinuities  in  waveguides,  the  quantity  which 
is  of  physical  importance  is  the  reflection  coefficient.  In  this  section  we 
calculate  this  for  the  fiindamental  mode  when  ka/n  <  1. 

The  reflection  coefficient  is  the  residue  of  ;2^(s,y)  at  s  -  -ik  and, 
from  Eq.  (7),  it  is  given  by 

0^-ik,d)  gg 

-nra —  *   ^   *  T  • 

We  have  from  T?q.    (n) 

«V-.  H^        T  1        hj-^       ^00     N.(-^X(^m>^) 

^^(-xk,d).  iim  ^j-^j_-j^*  j:  

s-o-ik     +1  ni=o        '^m  m 

and 


N   (-ik) 


N  (s)  '^  TTxIr      "  -it  a»  s  ^  -ik. 


Thus  the  reflection  coefficient  R  is  given  by 


Substitution  in  equation  (17)  for  p   (ik,d)   gives 
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—  *    L  (-ik)  ]  T'  ^ 


(16) 


V  d 


ig 


aR    K^^^^     j 

T'    L  (-ik)  ' 


aR 


-  1 


.  d 


For  aiir-plicity  it  is  now  convenient  to  introduce  the  non-dimensional  expressions 
H^(s),  H_(s)  defined  by 

If  this  change  of  notation  is  made  in  equation  (lb),  then,  writing  a  for  c/a, 
we  obtain 


(19) 


R 


Tl^^^"   a 


H^(ik)  ]   K^(ik) 


R  .  . 


1  - 


Hj(ik)  )   H^(ik) 


Equation  (19)  may  be  inverted  to  give 


{■^) 


R 


6(i+ig)H:(ik) 


i?(i 


a+H;(ik)-ig(a-r(ik)) 


Ig 


where  5  is  written  for  l-a»     The  expressions  in  Appendix  A  clearly  indicate 

iO 
that    |H   (ik)|   »  1  and  hence  that  H  (ik)  may  be  written  as  e        (with  ©  real). 

It  is  instructive  to  consider  the  physical  significance  of  the  expression 

defined  by  Eq.   (20)  when  g  -  0.     This  expression  may  be  derived  from  physical 

considerations,   as  will  now  be  shown.     We   assame  that  the   step  is  absent  and 

-ikx 
superimpose   a  field  Ae  incident  from  the  region  0<y<d,   x<0  on  the 

+ikx 
original  field.     This  incident  field  will  cause  a  wave   AR,e  to  be  re- 

flected back  in  the  negative  x-direction  (R,   may  be  calculated  by  use  of 
the  V;iener-Hopf  technique)   and  we  thus  find  that  the  total  velocity  component 
of  the  potential  in  the  fundamental  mode  for  y  <  d,  x  <  0  is 
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AILe     +  e     +  Ae 

and  we  choose  A  so  that  the  partial  derivative  of  this  expression  with  respect 

to  X  vanishes  on  the  line  x»0,0<y<d.  If  this  is  done,  then  the  trans- 

-ikx 
mission  coefficient  for  the  problem  of  a  wave  Ae     incident  in  the  region 

formerly  occupied  by  the  step  is  seen  to  agree  with  the  expression  of  Eq.  (20) 

for  g  =  0,  Thus  we  see  that  the  expression  (20)  illustrates  the  connection 

which  was  alluded  to  in  Section  I  between  the  bifurcation  and  the  step  problems, 

The  assumption  g  »  0  of  course  merely  implies  that  Y  is  constant  over  the  short 

vertical  part  of  the  step.  For  ka  <  n,  it  is  possible  to  compare  the  present 

results  with  those  given  in  the  handbook, and,  since  the  solution  given  there 

is  an  equivalent  circuit  formulation,  the  present  results  have  to  be  recast 

in  this  form.  The  equivalent  circuit  iss 

X 


^  1 "y'  y   , 

0  0  o    1 


xB  ~r 

^0 


where  Y  .  Y  are  the  characteristic  admittances  of  the  two  guides  and  B  is 
o'  o 

purely  real.  —  «  G  may  be  expressed  in  terms  of  the  reflection  coefficient 

y' 

o 
R,  and  we  have 

{(a+l)R+(a-l)} 
(21)  iC ^^^^ 

If  this  expression  for  R  in  terms  of  g  is  inserted  in  Eq.  (21)  then  we  obtain 

r22l        C  .  ^  J  ^i"Q  *   g  '^^^Q 
'■'^^''  ^  a    \   cose  -  g  sin» 

Inspection  of  this  expression  for  G  shows  that  it  maybe  obtained  by  applying 
the  above  simple  argument  for  determining  the  reflection  coefficient  to  the 
case  in  which  the  bifurcation  extends  a  distance  tan  g  to  the  right  of  the 
origin,  and  we  then  obtain  the  exact  form  of  C.  Thus  the  reflection  coeffi- 
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cient  for  the  step  problem  is  the  same  as  the  transmission  coefficient  for  the 
bifurcation  problem  mentioned  above,  when  the  bifurcation  extends  a  distance 
tan"  g  further  towards  the  positive  x-direction  than  the  original  step  stiructure , 


V.   Comparison  with  the  results  of  the  Handbook 

Some  simplification  wiDl  be  afforded  if  Eq.  (16)  is  rewritten  in  dimen- 

sionless  form  by  making  the  substitutions  given  in  Section  IV  and  defining 

-1/2 
a  new  set  of  constants  d  as  c  d    .  The  following  equation  is  then  satisfied 

n     n 

by  the  d  's: 
•^      n 

7r    00  d  (|i  li,  -6  ; 
(23)  H  (  K  )d  -  a  -  a  r  -"^  ^^      y     H  (-  \L) 

2    2  2 
where  e  =  kd/n  and  ij,  =  m  -e  .  The  function  g  may  now  be  written  in  dimen- 

1/2 
sionless  form  as  cfa    where 


m=l 


The  formulas  are  now  in  the  simpliest  form  possible  and  are  suitable  for  com- 
parison with  the  results  of  the  Handbook.  First  of  all  we  shall  consider  the 
special  case  of  neglecting  all  powers  of  e  >  2,  To  thj.s  order  of  approximation 
it  may  be  shown  that  for  m  >  0 

A"'ml  (m(l-6)/6)l 


«*(K„) 


where  1-6 

(v^)  ^ 


(ra/6)l 


The  argument  of  H  (ik)  may  be  shown  to  be  e  log  A,  and  we  thus  have 

^  =  ^  (log  A  +  a^/^f). 
e     a 

The  amount  of  computation  involved  in  this  special  case  is  fairly  small. 
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since  to  the  present  order  of  approximation  the  a  's  reduce  to  m.  Eqs.  (23) 
have  been  solved  on  the  assumption  that  all  the  coefficients,  apart  from  the 
first  five,  vanished.  This  gives   for  5  =  »,  the  result  C/c  (e  =  0)  =  1.556, 
while  the  corresponding  result  in  the  Handbook  is  1,36U.  For  6  =  .8  the  pre- 
sent  method  gives  . .■  -j  v  -   2.U5  while  the  value  in  the  Handbook  is  2.J48, 

In  order  to  test  further  the  accuracy  of  the  present  approach,  the  general 
expressions  developed  in  Appendix  A  were  used  to  calculate  H_(-  %S)   for  certain 
values  of  e,  and  an  approximate  solution  of  Eq.  (23)  was  obtained.  The  calcu- 
lations were  carried  out  for  6  =  U,  5  =  K-5e=  .6U,  6  »  .9;  and  the  values  ob- 
tained for  C  were  .95U  and  3.67.  The  values  given  in  the  Handbook  are  .96  and 
3.73.  The  values  of  e  and  6  were  chosen  so  that  ka/n  ranges  from  0  to  .3} 
hence  if  thei^  were  any  large  deviation  between  our  results  and  those  in  the 
Handbook,  it  would  have  been  evident  in  the  numerical  results.  It  is  thus 
evident  that  there  is  very  gcod  agreement  between  both  methods,  and  that  the 
solution  of  the  first  five  equations  is  sufficient  to  give  a  very  accurate 
answer.  In  fact,  complete  neglect  of  g  produces  a  percentage  error  of  less 
than  30  in  the  expression  for  C  while  the  percentage  error  in  R  will  of  course 
be  considerably  less,  (i.e.,  a  very  large  error  in  determining  g  does  not  pro- 
duce an  equally  large  error  in  the  expressions  for  R  and  C).  This  indicates 
that  the  assumrtion  of  a  constant  value  for  T  on  the  vertical  step  is  a  fairly 
satisfactory  one.  It  would  be  possible  to  obtain  an  estimate  for  the  error 
made  in  the  above  calculations  by  using  an  approach  similar  to  that  which 
will  be  indicated  later  for  e/6  >  1.  However,  since  the  above  numerical 
values  are  sufficient  to  indicate  the  accuracy  of  the  method,  it  does  not 
appear  worthwhile  to  enter  into  a  detailed  discussion  of  estimating  the  error 
for  this  case. 
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VI,  The  case  of  ka/n  >  1 

We  now  consider  the  modifications  which  have  to  be  made  to  the  above 

analysis  when  ka/n  >  1  but  kc/n,  kd/n  <  1.  One  difference  is  that  there 

will  be  a  term  in  T  which  behaves  s.3   exp(-i"nx),  where 

y.        ,.22/2  a/2 
y\  -  (  k  -n  /a  )  '  , 

as  X  -ooo  and  >}  is  now  real.     The  coefficient  of  exp(-i>»x)  ra^  be  fovind 
after  some  manipulation,  to  be, 

costy     2y^cv^     J,         2i(e-8)  ^    ^m^V^^M  ^  if2 

where       e     -dH(-K)^,         p-   (e^-e^)-""/^  and 


m        m  -         mm' 


(l-6)p 


91  .  -^S^(p,5,0)  +  S^(  iilii^  ,a,c)     -  S2[  §,1,0  j    -  S^(e,0,o: 

i^M  I  1  +  a  log  a  +  6  log  S 


and 


Sjj(x,a,p)  -  Yl    Y^^ 


00      ,  - 

-1  X  X 


y(n-?)^-c 


H^N  [  v/r7TT2~^       " 


The  S^'s  are  tabulated  in  the  Handbook, 

The  above  expressions  are  now  in  the  most  convenient  dimensionless  form  for 

numerical  computation. 

Another  different  feature  of  the  case  ka/n  >  1  is  that  the  functions 

H  C-  K   )   are  no  longer  real  and  we  have  H  (  K    )  =  T  (m-  +iP)G  (K„)  where 
—  ■       n  +       m         0       m     "^      +       m 

the  G    's  are  real.     The  equation  for  the  new  unknowns       e     defined  above  is 
+  n 


(2li)  ti  e  H  (K  )  =  a  -  c 

\   /         '^n  n  +   n 
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2,..  ^    I  ^     %^^m^n'^^^ 


^—       m     m  n 


Eq.  (2U)  may  be  written  as  two  sets  of  real  equations  as  follows 


(25) 


^^n  ,2,  ^  w  2  2.   ,,  4n^!^^n)  „  "?  ^  ^  ^^A'*^^^ 

-J-  G^(  K^)(^^-|3  )-  2p  -2 a  -  r.^m  (^  ^  u  ■) 

a5  05  m»l     '^m    n 


(26)     -r-  ( VP  ^^n  *  ^2 ^r^'-r^n     (u  .  u  ) 

a6  ao  in»»l      m    n 


where   e   «  f  +  ig  . 
n    n    ^n 

It  is  possible  to  reduce  Eq.  (25)  and  (26)  to  two  sets  of  linearly  in- 
dependent equations  but  since,  in  the  present  case,  they  were  solved  on  the 
digital  computer  there  does  not  appear  to  be  any  advantage  gained  by  their 
reduction.  The  above  equations  also  have  the  advantage  that  they  are  in  a 
form  from  which  it  is  fairly  simple  to  estimate  the  error  involved  in  their 
solution. 

VII.  Solution  of  the  equations 

In  this  section  we  shall  indicate  two  approaches  which  may  be  utilized 
to  estimate  the  error  made  in  solving  only  a  finite  number  of  Eqs.(25)  and 
(26).  It  is  of  course  practically  impossible  t o  make  a  precise  statement 
regarding  this  error,  due  to  the  complexity  of  the  equations.  Thus  some  of 
the  following  analysis  would  have  to  be  modified  in  certain  cases  but  the 
main  scheme  of  attack  will  always  remain  the  same.  Solving  only  a  finite 
number  of  equations  is  of  course  equivalent  to  assuming  that  the  coefficients 
on  the  right  hand  sides  of  EqE.(25)  and  (26)  are  non-zero  only  for  a  finite 
number  of  terms.  The  basis  of  the  approach  used  here  is  to  estimate  the 
error  made  on  the  right  hand  sides  of  the  equations  by  making  this  assumption 
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and  then  attempting  to  obtain  bounds  for  the  error  made  in  determining  the 
coefficients  by  using  the  known  estimate  of  the  error  made  on  the  right  hand 

sides  by  cutting  off  the  series  after  a  finite  number  of  terms, 

2     2  2 
It  is  clear  that  all  the  terms  of  the  form  (u,  a  -e  ),  (a  -6  )  occurring 

m  n       n  '^ 

2  2 
in  Eqs.(25)  and  (26)  will  be  positive  if  M--e  is  positive  and  thus  a  criterion 

for  all  these  terms  to  be  positive  is  that  e  <  1/  )^,     This  case  will  be  con- 
sidered first.  Examination  of  the  equations  clearly  indicates  that  it  is  not 
possible  to  make  an  exact  analytic  statement  concerning  the  forms  of  f  and  g 
and,  in  particular,  it  is  not  possible  to  assign  definite  signs  to  them.  It 
is  clear,  however,  that  they  will  ultimately  attain  a  fixed  sign  and  it  does 
not  appear  reasonable  for  them  to  change  sign  more  than  once  as  n  increases. 
In  any  case,  it  is  possible  to  say  what  combination  of  signs  of  f  and  g  can- 
not occur. 

If  it  is  assumed  that  the  f  's  are  all  negative,  then  Eq,  (25)  shows  that 
the  g  's  are  all  negative  and  then  from  Eq,  (26)  this  last  statement  implies 
that  the  f  's  are  all  positive.  Hence  the  initial  assumption  that  the  f  's 
are  all  negative  is  incorrect  since  it  leads  to  a  contradiction.  A  similar 
argument  shows  that  the  g  's  may  not  all  be  positive.  The  assumption  that 
the  f  's  are  positive  and  such  that  the  right  hand  side  of  Eq.  (25)  is  nega- 
tive leads  to  the  result  that  the  g_'s  are  all  positive  and  hence  a  solution 
with  the  f  's  satisfying  this  condition  is  not  permissible. 

The  predominance  of  the  left  hand  side  of  Eqs.(26)  and  (27)  suggests  that 
the  form  of  solution  most  likely  to  occur  in  practice  is  that  in  which  the  f_'s 
are  all  positive  while  the  g  's  sire  all  negative.  For  the  case  of  e  >  1/ V^ 

it  is  impossible  to  make  any  statements  whatsoever  regarding  the  signs  of  f 

2     2  2 

and  g     since  for  the   first  few  values  of  n  the  coefficients   (m-„-P    )*    (^^Jn^^n~^   ^ 

will  be  negative   and  the  above  argument  concerning  the  admissible  combinations 
of  signs  will  no  longer  be  valid.     For  the  purposes  of  the  subsequent  analysis 
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it  will  be  assumed  that  the  form  of  the  solution  is  such  that  all  the  g  's 

n 

are  negative  while  all  the  f   's,  apart  from  the  first  r,   are  positive.     The 
form  of  the  terms  on  the  left  hand  side  of  the  equations  suggests  that  the 
present  assiimption  is  the  most  likely  one  but  it  is  not  the  only  one  which 
may  occur.     It  has  the  advantage,   however,  of  being  sufficiently  simple  to 
illustrate  the  method  that  could  be  employed  in  any  particular  example  with- 
out making  the   analysis  too  elaborate  and  difficult  to  follow. 

The   above  arguments  break  down  in  certain  critical  cases  e.g.,  when  the 
right  hand  side  of  Eq.    (25)  becomes  negative  for  n  greater  than  a  certain 
value.     No  criteria  can  be  given  for  these  critical  cases  and  each  particular 
problem  has  to  be  investigated  separately.     However  the   general  approach 
given  here  may  be  modified  to  suit  any  particular  case  under  consideration. 
In  some  cases  the  work  involved  may  be  considerable.     We  now  proceed  with  a 
more  detailed  study  of  the  equations,     Eq,   (25)  gives  that,  foi    n  >  r, 

^j/^-P^)oliK^)  4       r 

(27)  C<f  ^       —    <a       -y     X     f 

^      '  n  72  -t-.      ran  m 

a5  m=l 

where  X       »   (n  u.  -e  )  /  (ti^+p.^), 
mn         ^m'^n  /      '^m  ^n' 

Defining  N  as  some  fixed  integer  >  r  we  have  for  n  >  N 

/pp\  1  ^  1 

JTOV?)"        n^(l-e^  /  N^)^/^(l-(p^e^)/N^) 

Hence  we  have 

(29)  0  <  f„  <  %  (a-^/'  -  t    W.) 

n  m=l 

where  C  is  a  constant  which  may  be  determined  from  inequality  (28)   and  the 

inequality  B13  for  G  (K    )•     Similarly  it  may  be  deduced  from  Eq,   (26)  that 

for  n  >  N 


(30) 


where 


-  19  - 


?pC 


rr 


^"+e 


N 


These  inequalities  for  f     and  g     may  now  be  us«d  to  estimate  the 


aeries 


CO 


00 


I^\n'm  '  ^^ 


^     g 


From  equation  (29)  we  hav« 

00 


■-N+1  ffl-N+1  m  Sfl. 


(31) 


<  C 


r       f  ^.^V^ 


T-2 


's^s 


(/x  +e  +ti^) 


()^?I7  *  Hg) 


) 


and  the  right  hand  side  of  this  equation  may  be  written  as 


where 


n 


-1/2  ^'n,      /Vnx     e^    In       ,,  -1  ^^il 
n  n     "-  -* 


'Sn 


ilg        Ucn5+N)*HsH^-6^  1     ti„        jv^TV^Tnij^ 


2     ^    -T 


log 


Unipnv^ 


log 


THiTTTv^ 


S  f  n 


2 


'nn 


f^^°^    (1^0(.;-"^)     *'7['"^^ 
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Similarly  finsin  equation  (30) 

00 


^  ^nm^  ^'\*^  fs^^s'^Sn 


where 


^^ 


CO 


h         2 


'1  - 


TT 


n^, 


n 


y  +  — r-  log 
n 


6  -ti^M^+nN 


(0, 


Sn      n^^-Hg 


2  -'^    -y  *  :?  ^"e  — R 1  ^°g  -r- 


-J-.-T 


S  y«  n 


CO 


(2e^.n^) 


nn 


n 


D^i; 


n 


n  (ti|,j+M.jj)         n  e 


*    ^n  Hi      .     -1  ^^■ 


?]  . 


It  is  now  convenient  to  rearrange  the  first  N  equations  of  each  of  Eqs,  (26) 
and  (27).  They  may  be  rewritten  as 
2..-  ^  r 


(32)      ^^ 


n  n  '^   n 


X  f  -  a 

mn  m 


-1/2 


00 


^r^  mn  m 

N+1 


G^(K  )  f    o  o        9  1   N  00 

(33)     ^  ^^^n^  VP^^n  *  ^^4'n\  *  f  ^mn^  '  '  ^  \m^  • 


If  the  right  hand  sides  of  equations  (32)  and  (33)  vere  known  exactly  then 
it  would  be  possible  to  obtain  the  exact  values  of  f^  and  g^.  Ifrifortunately 

this  Is  not  the  case  but  the  bounds  between  which  these  right  hand  sides  lie  are 

-1/2 
known  in  terms  of  f, •••f  .  For  equation  (32)  the  bounds  are  a  '  and 

1    r 

r 


,--L/2 


-  Jn  *  C  ^S^^S  ^Sn  * 
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while  the  corresponding  values  for  Eq.  (33)  are  0  and 

r 


''^'^' 


S^^S^Sn  • 


If  f  ,  g",  n  <  1  <  N  are  defined  as  the  solutions  of  Eq.  (32)  and  (33) 
when  the  series  on  the  right  hand  side  are  neglected,  then 

(iu)  ^"  J[i  (1^ -p  )(f  -f  )-2pn^(g,-g'')U  - r X  f-r^  (f -1 ) 

^  -.2       ^^n  "^n  ^         n     n       ^^n  ti  ^n   (  jr-.   rcn  m    4-     mn^  in     m 

a6         1  I  N+1  1 

G?(K    )    f         o     o  *  o  .1  <»  N 

(35)  -^<^^n(VP^(V^'n)^2P^n(V^;)|  '  "  I^^W  ^  ^.n^^-^)' 

If  e  and  V  are  defined  by  If  -f*l,  Ig  -g*|  then  Eq.  (3ii)  and  (35)  yield  the 
n     n  ''nn''nn' 

following  inequalities 

(36)       *  '". n  In -p^ie  <r  ix_iif„i+r  ix_i£  +  — "g:(kjv^ 

^   '        _2    n'  n  ^  '  n  —  tr-%  '  mn'  '  m'  4—  '  mn'  m    .^   +   n  n 
a6  N+1         1  a6 

G^(K  )     o  o      °o  N         2pn^   ^ 

(37)       -^   ■"  Li  lu  -B  Iv  <  5"  IX   llg  1+7"  IX   Iv  +  ^  g'^(K  )£   . 

^  ' '        ~2   ^n ' ^n  ^  '  n  -  tr-n  '  mn '  '^m '  4-  '  mn '  m    .2   + '  n'  n 
a5  N+1         1  06 

The  maximum  values  of  the  first  series  on  the  right  of  each  of  Eqs.(36)  and 

(37)  will  be  denoted  by  f  and  u  respectively.  Addition  of  these  equations 
then  gives 

f     ?   ?     ?  W!(Nn)  N 

(38)  (VV  YnK-^    \'<n-^-^^   '^n  *  -n  ^  f  '^mn^^m  ^  ^m^' 

If  the  coefficient  of  e  +  v  on  the  right  hand  side  of  Eq.  (38)  is  negative 

n        n 

then  the  inequality  is  trivial,  but  since  e^  and  v^  are  interchangeable  in  in- 
equalities  (36)  and  (37),  it  is   still  possible  to  obtain  a  non-trivial  form  of 
Eq.   (38),     If  the  coefficient  of  (e^  +  v^)  is  defined  as  p^^  then  application 
of  Eq.   (38)  to  itself  gives 
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(39)  '  -    '-  -    '   ""' 


p  (c  +  V )  <  f  +  co„+  5"  — - — ^T  +  CO  +  r"  |x  <,|(6„  +  v^)y . 

■^n     n         n     -      n         n    t-^     p       ]  m        m    %-,  '  mS'     S         S   ( 


If 


N     |X      I 


then  if  the  above  iterative  process  is  carried  out  an  infinite  number  of  times, 
a  definite  bound  is  obtained  for  (e     +  v  ).     Subtraction  of  Eqs.    (3U)  and  (35) 
yields  the  inequality 

^l^K^  •>     7  ^  I.  ^P^'n     2 

(UO)         "     ."  g  (LL^-B'^)(e  -V  )  <r     +  co  +T"   IX      |(6+v)*  ^G^(K)(e     +v). 

^      '  72  ^n  ^n  ^   ^  n     n  n  n  *~,  '  mn'     m    m  .2     +^   '■n       n         n' 

ao  ni=l  ao 

This  inequality  will  be  trivial  if  it  yields  a  result  of  the  form  (c  -v  )  <  (1+6) 

(t  +v   )  where  6  >  0.     However,  since  addition  of  Eqs.   (3U)   and   (35)  produces 
n     n 

the  above  inequality  with  the  coefficients  of  e  -v  ,   c  +v     reversed,  it  is 
^  •>  n     n'     n     n  ' 

always  possible  to  obtain  bounds  for  e  -v  ,   c  +v     and  hence  for  t     and  v„. 
•'     ^  nn'nn  nn 

If  r  ?<  0  then  the  bounds  for  e   ,   v    will  be  obtained  in  terms  of  f-,,.f     and 

n   n  1    r 

hence  equations  of  the  form 

r 

are  obtained.  This  new  set  of  equations  may  again  be  treated  iteratively.  It 
shoxild  be  noted  that  the  above  approach  is  considerably  simplified  if  ka/n  <  Ij 
in  this  case  there  is  only  one  set  of  equations  and  its  form  indicates  that 
the  \inknowns  are  positive. 

The  type  of  approach  indicated  above  may,  with  suitable  modifications,  be 
applied  to  estimate  the  error  in  any  particular  case  but,  since  some  of  the  in- 
equalities obtained  are  rather  crude,  the  following  method,  which  is  numerically 
more  complicated,  may  be  employed  to  give  a  better  estimate  of  the  error. 

The  system  of  Eqs.  (32)  and  (33)  may  be  written  in  matrix  notation  as 
A  X  »  y 


n  —  -^  n   4r-T  nS  S 
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where  x  and  7  are  respectively  the  transposes  of  the  row  vectors  (f, ...f^,g, o..g„), 
and        ^  ^ 

-•«•   00  -X-   C»  00  00 

and  A  is  the  matrix  of  the  coefficients  on  the  left  hand  sides  of  Fqs.  (32), 
(33).  If  the  elements  of  the  inverse  of  A  are  denoted  by  a ,  .  and  if  the  suffix 
i  denotes  the  i-th  component  of  either  coliunn  vector,  then 

2N 

Though  the  exact  values  of  y.  are  not  known,  their  upper  and  lower  bounds  are 
given  in  terms  of  f-,,,f  and  the  upper  and  lower  bo'onds  of  the  x.  may  be  ob- 
tained as  follows. 

If,  in  Eq,  (Ul),  the  maximum  values  of  y,  are  inserted  for  a .  .  negative 
and  the  minimum  values  for  a.  positive,  then  the  resulting  expression  gives 
a  lower  bound  for  the  x. 's.  Similarly  an  upper  bound  can  be  obtained.  Hence 
inequalities  of  the  form 

r  r 

y  b.  ,f .  +  c,  <  X.  <  d.  +  y   e.  .f . 

5^1  ^J  J   i  -  1  -  1  4^^  ia  j 

are  obtained  where  the  b's,  c's,  d's  and  e's  are  all  known  in  terms  of  f, ,.»f  • 
Either  an  iterative  or  some  similar  process  may  now  be  used  to  determine  upper 
and  lower  bounds  for  f, ,..f  .  These  bounds  may  then  be  used  to  obtain  bounds 
for  the  remaining  f's  and  g's. 


VIII,  Numerical  results 

The  set  of  equations  occurring  for  the  case  ka/n  >  1  was  solved  by  assuming 
that  the  inteper  N,  referred  to  in  the  previous  section,  was  equal  to  5.  Since 
there  is  no  particular  physical  significance  in  the  actual  numerical  values  of 
the  coefficients  f  and  g,  ,  they  have  not  been  tabulated.  One  quantity  which 
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is  of  particular  interest  is  the  reflected  energy  crossing  a  plane  normal 
to  the  X  axis  which  will  be  proportional  to 


r." 


liiri   I  T*gdy, 


where  the  asterik  indicates  the  complex  conjugate.*  If  the  coefficient  of 

"■iv)  X 

e       I    ,  referred  to  in  Section  VI,  is  written  as  ILcosny/a,  then  the  above 


integral  may  be  evaluated  to  give 


-la  i 


k|R|2  *  1     >^  |R^|2 


Hence  the  (reflected  energy/incident  energy)  ratio  is 

2     2    3      2       2 
The  two  quantities  |R|  ,  |r|  +  |t  I%I  '   \'^2^   »  ^^   given  in  Table  I.  The 

numerical  values  show  that  the  ratio  of  the  reflected  energy  to  the  incident 
energy  increases,  for  constant  step  height,  as  the  waveguide  height  decreases. 
This,  of  course,  was  to  be  expected  on  physical  grounds.  The  percentage  of 
the  reflected  energy  which  is  due  to  the  |R-  |  tends  to  increase  as  the  wave- 
guide height  decreases  for  constant  step  height.  This  percentage  is,  however, 
fairly  constant  ranging  between  19 /o  and  2$/o  for  all  the  values  given  in  the 

table. 

2 

The  values  obtained  for  |R|  were  compared  with  those  obtained  by  assuming 

that  all  the  coefficients  of  the  set  of  equations  vanish,  and  the  two  results 

were  found  to  differ  by  only  a  few  percent  in  all  cases.  The  difference  be- 

2 
tween  the  two  values  obtained  for  |R- |  was  greater,  but  this  contributes  only 

a  small  amount  to  the  total  energy.  It  thus  appears  that  a  very  good  estimate 

may  be  obtained  for  the  reflected  energy  by  assuming  a  constant  value  of  the 

velocity  potential  on  the  step. 

The  general  methods  of  the  previous  section  were  not  used  to  estimate  the 
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errors  in  all  the  sets  of  equations  concernedo  However,  in  all  cases  the 
values  obtained  foi-  f^  and  g^  did  not  differ  by  more  than  5/o  from  those 
obtained  by  assuming  that  f-  and  g^  were  the  only  non- vanishing  coefficients, 
For  the  case  of  e  ■  ,9  the  general  methods  were  used  and  the  results  were 
shown  to  oe  correct  to  within  about  5/o. 


c  »   . 

9 

e  «   .8 

- 

e  -   .7 

6  -   .6            1 

ka/n 

1.8 

1.6 

I.U 

1.7 

1.6 

I.U 

1.2 

1.6 

I.U 

1.2 

1.1 

1.5 

I.U 

iKl^ 

.25 

.32 

.U 

.23 

.26 

.3U 

.U7 

.20 

.28 

.39 

.U9 

.17 

.20 

H.l^ 

.31 

.Uo 

.52 

.29 

.33 

.UU 

.58 

.25 

.35 

.31 

.63 

.21 

.26 
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Appendix  A 

In  this  appendix  formulas  are  developed  which  may  be  used  for  computing 
the  functions  N  (K  )•  li"  the  same  method  is  employed  as  that  used  by  Carlson 


and  HeinsL  J  then  it  is  fairly  easily  shown  that 


XT  f.^   -   1   ./^        P(6)    ,^  sjaloga  -  clogc  -  d  logdj 


±-o 


N  (-s)N  (s)  »  -1  where 


p(6)  -TT  i 

n«=l 


'      2  \7 


(i-^fr)   ^^n\  ^^[-Ife] 


It  is  not  really  necessary  to  retain  the  exponentials  in  the  above  in- 
finite products  when  the  complete  expression  for  N  (s)  is  being  considered, 
but  they  must  be  retained  for  evaluating  each  separate  product.  For  generality 
we  shall  consider  the  evaluation  of  each  of  the  products  separately  since  pro- 
ducts of  a  similar  nature  occur  in  the  other  problems  of  diffraction  theory  and 
only  slight  modification  of  the  following  analysis  will  be  reqxiired  to  cover 
other  cases. 

It  is  clear  from  the  expressions  occurring  in  the  text  that  we  must 
compute  expressions  of  the  form 


Though  most  of  the  computations  involved  in  this  work  were  carried  out  on  a 
digital  computer  it  is  also  possible  to  obtain  fairly  convergent  expansions 
in  powers  of  e  for  the  above  products.  These  expansions  are  convergent  for 
all  the  powers  of  e  and  6  that  occur  and  have  the  advantage  of  being  very  much 
more  suitable  for  numerical  computation  than  straightforward  multiplication  of 
the  terms  of  the  product.  In  order  to  obtain  the  numerical  values  for  a  large 
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number  of  values  of  e  and  6,  the  digital  computer  has  a  distinct  advantage 
over  any  other  form  of  computer.  The  program  developed  for  the  present  prob- 
lem has  also  been  used  for  computing  infinite  products  for  which  the  approach 
used  here  may  not  be  used  (in  these  products  the  nr  in  the  infinite  product 
was  replaced  by  the  n-th  zero  of  a  Bessel  fvmction  or  of  a  transcendental 
equation  involving  Bessei  and  Hankei  functions).  The  following  formulas  are 
suitable  when  not  too  meiny  values  of  e  and  6  are  required,  and  the  computations 
for  ka/Tt  <  1  were  carried  out  using  these  formulas.  Expansion  of 

y    n  Ti 
2 

the  following  result  for  m  j<  0 


and      n     in  powers  of  e   ,  together  with  some  algebraic  manipulation,  gives 


if-ft^V-'fj-E^] 


where  xl  is  the  factorial  function,   j  is  Euler's  constant  and 

P^(5)   -  a^e^  +  a^e^  *  a^t^  +  a^6^°  *  a^e^^ 
where 

8m6^a^  -  <(3),  USm^S^a^  -  -3m\(5)   -  25^^(3) 

2^6"^ a^  -  -5ra^(7)  -  3ra^5^^(5)   -  36^(3) 

5x2V5^a^  =  -35m^^(9)  -  lljOm^5^^(7)   +  36in^6'^(5)   -  205^^(3) 

3x2^V6^^a3-  -  63m\(ll)-35mV^.(9)-31m^5\(7)-30m^6\(5)-356\(3) 

where  ^(s)  is  the  Riemann  Zeta  function  i.e., 

00    ^ 
n-1  n 
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The  other  two  products  in  N^(  \  )  are  of  the  same  form  as  the  above 
with  a  replaced  by  c  and  d  respectively  and  hence  the  above  formulas  can  be 
extended  to  cover  the  other  two  products  by  replacing  6  in  the  expressions 
for  a-»»««aj^  by  y-^,  and  1  respectively.  These  results  are  valid  for  r  <  1 
but  may  be  extended  to  include  the  case  of  1  <  t?.  <  2,  e  <  1  as  follows. 

For  T  >  1  the  first  term  in  the  above  infinite  product  will  be  complex 
5 

and  hence  if  we  consider  the  integer  n  as  now  ranging  from  2  to  infinity 
this  new  product  may  be  treated  exactly  as  before »  The  ^(s)  in  the  formulas 
for  a^,.,aj-  should  now  be  replaced  by  J^(s)-1. 

The  case  of  m  =  0  requires  special  consideration.  In  this  case  the 
products  occurring  in  N  (ik)  may  be  written  in  the  fonn 


exp  i  I  tlf^in'^  H?  -  iJj 


00  J. 


The  functions  occurring  in  this  exponent  have  been  extensively  tabulated  in 
the  Handbook. 

Appendix  B 

In  this  appendix  we  obtain  some  inequalities  for  the  G  ( J(  )  of  Section 

+   m 

VI.  Gjlf.)   is  given  by 

G  (  «  )-  exp  -l-^  log  ^--  ♦  ^  log  -  -  T 

where  /  1     \ 

Each  of  the  above  infinite  products  satisfies  the  following  type  of  in- 
equality 


>^m 


-  29  - 


1 
2  .? 


(B2) 


n«s  I       5  n  J  '-      -l       n=s  I  6  n  I  L     '  J 

,  ^  L   "« 1    r  ""i 


and  hence  G  ( )(    )   satisfies  the  following  inequality 
u   (1-5) 


rn PT 0 


■)I 


^^ '^ J¥  ^°s  T§5  *  ^  l°4     <VKJ 


T^ 


(B3) 


ti„(l-5) 


m  0 


-)  I       exp  n 


m 


1-6  ,  6         1  1        1 

-r  ^°^  TTe  *  5  ^°s  6 


&l(l*  !^) 


T^ 


m(1)(i)    m1^^  i') 

in  m         1-0 


where 


C'(^)    -    ^.  C^  -ffiT^S^ 


00 

n=r  "^m 

In  order  to  discuss  Eq.   (B3)   it  is  necessary  to  establish  inequalities 
for  the   factorial  functions  and  the  infinite  products  which  occur  therein. 
We   shall  first  of  all  consider  the  factorial  functions,   and,  for  brevity, 
shall  define  a  new  function  f(|i  ,6)  by 


nr 


"^^ '^  ^'■^  1^  *  i  ^°e  \\     • 


If  we  apply  Binet's  second  formula    (5J    for  the  logarithm  of  the  factorial 
function  to  the  factorials  of  Eq.    (BU)  then  f(^j^»6)  may  be  rewritten  as 


f(li^,6)   -  exp 


(B5) 


log  2m 


,Jl.6).j(m).l*(|!^-5—  *  ^)log(l  *  ^;Jj=St) 


-  (I^  .  |)log(l  .  i-)  *   (V  7)1°^(1  *  ^) 


exi:h(in,6) 


where 


and         0  < 
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a  (1-6)  ,  u 


B, 


B, 


1  P  "^ 

53F    -    ^  <  X(x)  <    23?  for  X  >  0 


B     «  i,  B_  *  irjY     are  Bernoulli's  numbers. 
If  the  inequality 


X  - 


x^  2         3 

-J-     <  log(l+x)   <  X  -  ^  +  ^ 


is  applied  to  the  logarithmic  terms  of  B5  then  the  following  inequality  may 
be  established  for  h(m,6) 


(96) 


7       T?       77/ T  ^  iUi"*<J/   • 

■     ^     ■»•"§ 

I  ^T^j^y 

,   1     f       6^          ,    6^   .    1 

1 

6n^ 

< 

1     6^ 

^U-(-)^        "         ^^J 

(1-6)' 

v.. 

♦  1' 


Similarly,  application  of  the  inequality  for   ^(x)  gives  the  following  inequality 
for  j(m). 


^1/5       ,,    .A       ^     /    5' 


^  (j§^  .(1-6))-  -^  (-ii^  *  l).  !^(l  .  -£-^<  J(m)<  ^^  .  (1-6)) 
^^'in     -^^  ^     12nf  ^(1-5)^        ^     12^^^         (1-6)^  ^V^"^  ' 


m 


m 


(37) 


a,5^     Bpe^     Bpe^ 


2^;         12n;;       kil 
m  mm 


A  combination  of  Eqs.   (B6),   (B?)  now  gives  the  following  inequality  for  h(ra,5) 
(B8)  ^  log  2Ti|x^(l-e)+  h"(m,6)  <  h(m,5)  <  h''(m,6)  +  ^  log  2nHj^(l-6) 

where  h"(m,5)   is  obtained  by  addition  of  the  left-hand  sides  of  Eqs.    (B6),    (B7) 

and  h   (m,5)   is  given  by  the  sum  of  the  right-hand  sides  of  these  equations. 

In  order  to  complete  the  discussion  of  equation  B3  it  is  necessary  to  estao- 

lish  a  lower  bound  for  the  infinite  products  occurring  there.     They  are   all  of 

the  form 
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(B9) 


n"S  V      ^^m    '      n»s     ^      '^m 


The  following  inequality  holds  for  the  series  of  B9 


'n*s 


m 


m 


.00  2 


The  vsuLues  of  s  that  occur  are  either  1  or  2j  if  s  ■  2,  then  the  above  in- 
equality is  valid  if  s  is  replaced  by  unity  in  the  integral  and  if  the  first 
term  on  the  right  hand  side  is  omitted.  We  thus  require  an  inequality  for 

00     ^     ^"^m 

It  is  easily  verified  that,  after  a  little  manipulation,  the  following  result 
is  valid 


t'°<'-^4-^)'^- 


o  2  2 

"'00  ^^m    ' 


m  +  6  +  3 


Hm 


6-5 


e 

Tm 


The  right  hand  side  of  this  inequality  will  be  denoted  by  -c(ni*5)»  Finally, 
combination  of  B8  and  Bll  gives 


(1  -  L)   ^  r'^^f"^)  exp  J  h-(m,5)-  1^(»,6) 


^    <G>m) 


(B12) 


<  ^  V^^;;^^-^^     exp  J   h^m,5)  *i(m,l)   *  ^jt{m,   g|5:)-log(l  -  t-^Tt) 


-  log  (1  - 


e^(l-6) 


We  shall  only  require  these  inequalities  for  m  >  N  and  it  is  thus  possible  to 
obtain  simpler  forms  for  this  inequality.     In  the  exponent  on  the  right-hand 
side  of  Bi2  all  the  terms  are  positive  and  hence  the  inequality  is  still  valid 
if  all  the  terms  in  the  exponent  are  replaced  by  their  values  for  m  «  N;   also 
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2 

u  may  be  replaced  by  in(l  -  ^y)  '  •  In  the  exponent  on  the  left  hand  side 

of  B12  all  the  negative  terms  may  be  replaced  by  their  values  at  m  =  N,  while 
the  positive  terms  may  be  replaced  by  zero,  a     may  now  be  replaced  by  m. 
Thus  it  is  always  possible  to  establish  a  fairly  simple  ineqiSLlity  of  the 
form 

(B13)  ^  <  ^(^J  <  V        ""  -  ^ 

m  '  m  ' 

where  A  and  B  are  constants  independent  of  m. 
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